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Abstract 

An unusual photo-galvanic effect is predicted on the topological insulator surface when its semi- 
metallic electronic spectrum is modified by an adjacent ferromagnet. The effect is correlated 
with light absorption in a wide frequency range (from a few to hundreds of meV) and produces 
a pronounced response that is not only resonant to the photon energy but also tunable by an 
external electrical bias. The exceptionally strong peak photocurrent of the order of /xA/cm may 
be achieved at elevated temperatures with the illumination power of 1 W/cm 2 in the THz range 
on Bi2Se3. These advantages could enable room-temperature detection of far- infrared radiation. 

PACS numbers: 72.40,+w, 72.25.Fc, 78.68.+m, 72.80.Tm 
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To date the photo-galvanic effect (PGE), i.e., direct current generation via light absorp- 
tion in a structure with broken inversion symmetry, has been used as a unique experimental 
tool to study kinetic phenomena. Three different mechanisms are known to trigger the PGE: 
(i) the drag of electron gas by photons,-! 2 - (ii) the imbalance in the excitation of electrons 
with opposite momenta,- 1 ^ and (iii) the shear of electron non-equilibrium distribution in the 
momentum space in the course of relaxation to equilibrium.-^ In this regard, the spin-orbit 
interaction leads to a specific manifestation of the latter two cases as it has been demon- 
strated previously that the transfer of optical spin polarization can realize a direct current .- 
However, the induced current is generally of low density due to the secondary role of spin- 
orbit coupling in the formation of band structure. Other approaches of PGE invariably 
suffer from the same fate of feeble response. 

One potential exception may be the topological insulators (TIs) — the material with 
spin-orbit interaction so strong that it inverses the conduction and valence bands and ex- 
hibits numerous unusual phenomena deviating from the conventional perception about the 
magneto-electric effects.- - — It is plausible to expect the peculiar properties of TI surface 
states also modify the PGE, particularly when the time-reversal symmetry is broken by 
external fields. For instance, absorption of circularly polarized light was recently shown to 
generate a direct current in the surficial region of the TI.— Generally, this circular PGE 
is conditioned by the spin-momentum interlocked nature of the TI surface states. To get 
a nonzero response, however, a magnetic field B must break the three-fold rotational sym- 
metry that results in an extremely weak photocurrent of ~ 10 fA/cm under the excitation 
power of 1 W/cm 2 and B = 10 T. A more pronounced effect (~ 10 nA/cm at T = 15 K) 



was reported in Ref. [12J, where the TI surface isotropy is broken by the obliquely incident 
light with an excitation energy much higher than the TI band gap. While the spin-orbit 
interaction is the key to the circular PGE, the weak low-symmetry interactions limit the 
density of dc transport to a low level in both cases as before. 

In this study we draw attention to another, much stronger mechanism of PGE that 
stems from the band structure asymmetry at the interface between a TI and a proximate 
insulating or dielectric ferromagnet (FM). The basic principle of this effect is illustrated in 
Fig. 1. Here, the FM is assumed to be magnetized along the y axis that results in the shift 
of the TI surface Dirac cones along the x direction away from the V point in the momentum 
(i.e., k) space.- 1 ^ We focus, however, on one consequence of this shift beyond the contact 
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point that is manifested in the contortion of the conic shapes due to the spin-independent, 
quadratic-in-A; term [see the dashed curves in Figs. 1(b) and 1(c)]. Specifically, the states 
on the right side of the shifted Dirac cones in the figure are farther from the T point, 
in which the k 2 term in the Hamiltonian increases more rapidly for both conduction and 
valence bands. For the n-type TI plotted in Fig. 1(b), the optical excitation is suppressed 
for the states on the left side of the Dirac cone, since both initial and final states are 
situated below the electrochemical potential \i. On the other hand, the right side of the 
cone allows the transitions between resonant states via the same photons because they (i.e., 
the electronic states) are located astride \x. Consequently, a significant imbalance induced in 
the momentum (i.e., group velocity) distribution of the excited electron ensemble produces 
a net photocurrent. Similarly, the p-type TI [Fig. 1(c)] gives a nonzero photocurrent along 
the opposite direction because the optical excitations on the right of the Dirac cones are 
suppressed instead. 

As a didactic but practically important example of the approach, we start with a discus- 
sion on optical absorption in the TI/FM heterostucture before proceeding to the PGE. The 
corresponding states are described by the effective Hamiltonian of the low-energy TI sur- 
face electrons in an exchange interaction with a proximate magnet. Keeping the linear and 
quadratic terms of the two-dimensional electron momentum k = (k x , k y , 0), the Hamiltonian 
takes the form^ 



where z = (0, 0, 1) is the unit vector normal to the TI surface, vf is the Fermi velocity, a = 
(a x , (T y ,a z ) is the vector Pauli matrices of the electron spin, and D is the material parameter 
for the quadratic term. The last term in Eq. ([I]) describes the exchange interaction with the 
FM in terms of the effective exchange field G = £M, where M denotes the FM magnetization 
and £ is proportional to the exchange integral.— 1 ^ Since the exchange energy involves only 
the electron spin operators, its effect can be simply regarded as a momentum renormalization 
for the first term in Eq. (TO , with k y — » k y + G x /Uvf and k x — > k x — Gyjhvp. In the case of 
M directed along the y axis (i.e., G y — G), the eigen-energies of the Hamiltonian are 



where p =Hvf^- is the momentum (in energy representation) of the conduction (b = c) and 
valence {b = v) electrons and d = D/h 2 v F . Note that the contact point energy corresponds 



H = hv F [a x k] • z + Dk 2 



(1) 




(2) 
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to e CtP = e V)P = dG 2 due to the quadratic term. Accordingly, the relations e CtP > dG 2 and 
£u,p < dG 2 are satisfied by the conduction and valence bands, respectively. The reference 
(i.e., zero) of energy is defied at the Dirac point of unmodified TI surface states. 

Let us consider a normally incident electro-magnetic wave whose electrical compo- 
nent E generates the electric dipolar transitions. The case of monochromatic wave 
E = Eoeexp [i(q ■ r — cut)} corresponds to a certain polarization, where u is the angular 
frequency and q = qz is the photon wave vector. The circularly polarized wave is described 
by polarization vector e = e± = (x ± iy)/y/2, while e# = (xcos# + ysin^) corresponds to 
linear polarization with the angle 9 between the polarization plane and the x axis. The effect 
of unpolarized light can be described as a half sum of the contributions induced by non- 
coherent waves with perpendicular linear polarizations or opposite circular polarizations. 

In the Coulomb gauge, the vector potential for the electro-magnetic radiation has the 
form A =Re(icE/w), where c is the velocity of light. From the perturbation treatment, the 
Hamiltonian of electron interaction with the radiation is 

evFr^ . , „ 2eD . 
V = — - [a x A] ■ z + — — k • A , (3) 

c nc 

where e is the electron charge. Neglecting the small term containing D, the squared absolute 
value of electron-photon matrix element is given as \M\ 2 = v 2 7 e 2 E' 2 /&uj 2 for both e + and 
e polarization.— Thus, this is also applicable to unpolarized light. The effect of linear 
polarization appears in \M\ 2 as the additional factor 2 sin 2 (<^ p — 6) (ip p is the angle between 
p and x), which reduces to 1 after averaging over 9. An important property of the matrix 
element M is its independence on the electron momentum and exchange interaction. 

The rate of photon energy absorption P = hw^2 (W VtP - CtP — W CtP - v , p ) is associated with 
the probability of interband transition \b, p) — > \b',p): 

W b ,p;v,p = — |M| 2 /o(^, P ,/u)[l - fo{ev,p,fj)]S{\e VtP -e btP \ - hu). (4) 

Here we ignore the small photon momentum q compared to k as well as a minor deviation of 
the population factor f(v, p) from the equilibrium Fermi-Dirac function fo{£v,p, A*)- Applying 
Eq. (J2J) to Eq. (j4j), straightforward integration of H^, jP;CjP — W c>p;v>p over p expresses the 



absorption rate P in terms of the roots of the energy 5 function p x = p x = G± « / (hu/ 2) 2 — p 2 
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(A = ±) as well as the residual integral over p y : 



hw/2 



P(/^,/i) = |p f E 



/o(& 



'«,PA> ~~ /o (^c,p A , ^) 



dp. 



(5) 



-hui/2 



where a = e 2 /hc is the fine-structure constant and Pj = (c/8tt)Eq is the energy flow of the 
incident light. This equation allows simplification under the assumption of relatively small 
energy shift dG 2 compared to the thermal energy ksT. The final output in the normalized 
form (i.e., normalized absorption) can be written as 



determines the absorption ratio. At low temperature [y » 1) in the absence of band cur- 
vature (d — 0), this value is four times smaller than the graphene absorbance tra ~ 2.3%,— 
reflecting the absence of degeneracy on valleys and spins. It is also obvious that 7ra/4 is the 
maximal structure absorbance attainable with a small \fj,\ and high photon energy- 
Extending the analysis to the PGE begins with the general expression for the surface 
current in the diffusive regime 



where the axis x is directed along the electron flow, Aq is the area of the TI/FM interface, 
and v(b, p) = v FV p £fc iP is the electron group velocity. The electron distribution function 
f(b,p) contains the non-equilibrium effect due to light absorption [i.e., Eq. (jlj)]. When 
f(b,p) = /o(A,p; /-O (i- e -> in equilibrium), the current j becomes zero exactly. This also 
means the absence of dark current independent of the energy band distortions. Accordingly, 
we hereinafter focus on the deviation from the equilibrium distribution; i.e., A/(6, p) = 
f(b,p) - /o(e6, P> ^)- 

In the stationary case, the relaxation time approximation can yield a simple solution to 
the kinetic equation for the distribution function as A/(6, p) = S(b, p)r(6, p), where r(b, p) 
is the relaxation time and S(b, p) is the rate of electron generation in the state \b, p). In turn, 
the generation rate can be expressed via the probability of transitions discussed earlier [see 
Eq. (J4])]; i.e., S(b,p) = J2v (W&',p;&,p — Wf>,p;&',p)- To extract the explicit effect of excitation 
imbalance in the electrons with opposite group velocities, it is sufficiently to restrict the 



Pi 




in terms of arguments v 



hu/2kBT and n = (4/x — dh 2 uj 2 ) /AksT . The factor ^7raF(is,ri) 




(7) 
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consideration to a constant relaxation time r(6, p) = r independent of momentum p and 
energy band b.— As such, the population imbalance is solely defined by photo-transition rates 
S(b,p). At a fixed p y and the photon energy Jtw (> 2 |p y | ) , the only states with p x = p± 
are involved in the interband transitions, while the proximity effect produces the energy 



difference £b,p+ ~ £ &,p- = 4<iGy (hu/2) 2 — p 2 for both electron and valence bands (Fig. 1). 
This asymmetry between £b,p + and £f> iP _ presents the main reason for direct photocurrent 
generation. Indeed, it can be shown after some algebra that the summation over p x in Eq. ([7]) 
produces an expression determined by the difference of the population factor rather than 
the velocities of photo-generated electrons and holes. The net photocurrent effect gives the 
integral over p y for a given hoj provided that |A/(6, p)| <C /o(£&, P 

frw/2 

J = 2 1x K i v F ^ M ^ J dPy ^ ^° (^' p+,/i ) ~ f° ( £ «.p-»^)] ~~ [fo ( £ c, P+ ^) - fo ( £ C ,p_,/i)] } • 

-hw/2 

(8) 

Along with numerical evaluation, an analytical approximation of j would be desirable 
for a qualitative analysis of the effect. Assuming that the temperature is not too low 
(ksT > dG 2 ), Eq. (J5]) can be represented in terms of the incident light power and function 
F: 

J^-™ dG ^ P i F MF(r,,v), (9) 

which explicitly displays the dependence of the PGE on absorption [Eq. (jBJ)]. This ex- 
pression clearly illustrates that the photocurrent originates from both the proximity effect 
(~ G), which breaks the isotropy of the interface Hamiltonian, and the spin independent 
parabolicity (~ d) of the energy bands. Furthermore, the dependence on photon energy 
(i.e., dimensionless variable v) manifests a well-defined maximum that can be tuned by a 
change in the chemical potential. Note that switching in the FM magnetization (G — > — G) 
leads to reversal of the photocurrent, which is in accordance with time-reversal symmetry 
inherent in the linear PGE.- 

The photocurrent calculated in the structure with a typical TI is shown in Fig. 2 for 
different excitation energies at 300 K and 4 K. Specifically, Bi 2 Se3 is assumed along with 
t = 1 ps, G = 40 meV, and the incident radiation power density of Pi = 1 W/cm 2 . 
The maximum current density at room temperature is of the order of 0.1 /iA/cm, whereas 
it reaches over 10 /zA/cm at liquid helium temperature. These numbers are drastically 
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larger than the photocurrent reported earlier in TIs.— ^ Another remarkable property of 
the proposed PGE is its correlation with the normalized absorption curves. As shown in 
Figs. 2(a) and 2(b), stronger and broader absorption characteristics leads to a similar trend in 
the photocurrent response at room temperature. On the other hand, the sharp transitions 
in the absorption curves (vs. fi) at low temperatures lead to the corresponding resonant 
current peaks [see Figs. 2(c) and 2(d)]. This can be understood clearly from Figs. 1(b) and 
1(c); note also the relation dP/drj oc j from Eqs. ([6]) and (Q although its validity is limited 
at low temperatures. 

As indicated earlier, the main contribution to the photocurrent comes from the asymmet- 
rical distribution of the photo-excited carriers that is in turn due to the shift G of the Dirac 
cones and the Dk 2 term in the Hamiltonian. While the parabolicity is proportional to the 
parameter D, which is well-known for most TIs, the magnitude of the exchange energy may 
vary depending on the TI/FM materials as well as the quality of the interface. Based on the 
earlier studies in the literature,-^— it appears that G = 5 — 50 meV can be expected. In 
the limit of low G, Eq. (J9j) predicts proportionality of the PGE to G. However, numerical 
evaluation of the more general expression [Eq. ([8])] illustrates a saturating trend in the peak 
magnitudes of resonant curves j(fi) as G increases (dG 2 3> &bT). This is clearly visible in 
Fig. 3(a), where increasing G from 5 meV to 40 meV at T = 4 K significantly broadens the 
resonant photo-response but introduces only a minor change in the peak height. The satu- 
rating behavior can be explained, at least in part, by the corresponding photon absorption 
curves. As can be seen from Fig. 3(b), an increase in G leads to a more gradual change in 
P with a decreased slope. Since both factors (i.e., G and dP/drj) contribute, their impacts 
compensate each other and the peak current remains largely unaffected. 

Figure 3(c) shows the position (in /i) of the peak current as a function of the photon 
energy. The observed dependence is nearly linear in a wide range of photon energy, suggest- 
ing the equivalence (or strong correlation) between the frequency scanning and the chemical 
potential (thus, gate bias) sweeping in the detection of the PGE. In other words, the energy 
of incoming photons can be determined by measuring the TI surface current as a function 
of the gate bias. The strength of the photo-response is particularly pronounced at low tem- 
peratures as shown in Fig. 3(d). Note that all curves are calculated assuming an identical 
relaxation time r = 1 ps. As the relaxation rate tends to decrease rapidly with temperature, 
further enhancement of PGE beyond tens of /xA-cm/W appears likely at T = 4 K. Moreover, 
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the relaxation time can be much longer than 1 ps even at room temperature as the sample 
quality improves. 

The strong PGE discussed above can be used to realize a photo-detector in the THz/far- 
infrared frequencies. One clear advantage of this mechanism is that the dark current can 
be eliminated without sacrificing the response time since the resulting photocurrent does 
not require an external bias. In conventional photo-detectors (such as those based on p-n 
junction photodiodes, Schottky barriers, quantum well structures, etc.), an external bias is 
applied to accelerate the excited carriers or to decrease the pn junction capacitance for fast 
response. However, the consequential dark current is detrimental to the device performance 
at room temperature as it is a powerful source of noise. The prospective device, on the 
other hand, is expected to work at room temperature with an extraordinary sensitivity, 
while having a very short response time. Since no junction is formed in the channel, the 
response time is determined by the carrier transit time. Not only the excited carriers have 
a high initial ensemble velocity (~ vf) similar to the Schottky-barrier detector, but also the 
momentum relaxation in ideal TI surface states is expected to be very long, particularly for 
low energy electrons, due to the suppression of backscattering. Overall, the structure offers 
a unique opportunity for room-temperature detection of long-wavelength photons. 

In summary, a linear PGE on the topological insulator (TI) surface adjacent to an insulat- 
ing FM is examined. The dc photo-excitation stems from the TI surface band modification 
under the influence of both the symmetry-breaking proximate exchange interaction and spin- 
independent parabolicity of the dispersion law. The phenomenon is invariant to the sign of 
light circular polarization, thus effective at unpolarized excitation. At the same time, the 
PGE reveals the frequency-dependent resonance-like features, particularly for low temper- 
atures, despite the absence of any discrete energy levels. The strong correlation between 
the excitation energy and photo-response enables straightforward detection of radiation fre- 
quency via electrical means such as the gate bias.— The photocurrent estimated with a 
relatively short relaxation time of 1 ps exceeds those of the previously reported PGEs by 
orders of magnitude. The proposed TI/FM structures may have significant advantages over 
the conventional devices in the detection of long-wavelength photons beyond the thermal 
noise limit.— ^ 
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(a) 



y y " 




FIG. 1: (Color online) Schematic illustration of photocurrent generation on a TI surface adjacent 
to an insulating ferromagnet (FM). (a) TI/FM hybrid device structure. The arrows indicate the 
incident electro-magnetic radiation and the block arrow denotes the magnetization of the FM layer. 
(b,c) Photocurrent generation in a n-type TI and a p-type TI, respectively. The dashed lines are 
the TI surface bands taking into account the k 2 term and the proximity effect, while the solid lines 
only include the proximity effect. The shaded area denotes the filled states under electrochemical 
potential \x. 
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FIG. 2: (Color online) Photocurrent density and normalized absorption curves at (a,b) 300 K and 
(c,d) 4 K versus the chemical potential with different excitation energies:— hcu = 4.4 meV, 13.7 
meV, 50 meV, and 135 meV for curves 1, 2, 3 and 4, respectively. All calculations assume the 
relaxation time r of 1 ps and the proximate exchange energy G of 40 meV along with the material 
parameters corresponding to Bi2Se3 {vf = 4.28 x 10 cm/s and D = 13 eV-A 2 ). 
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FIG. 3: (Color online) (a) Resonance-like photocurrent and (b) normalized absorption curves at 
4 K with the excitation energy of 135 meV. Curves 1, 2, 3, and 4 represent the results with the 
proximate exchange energy G of 5 meV, 10 meV, 20 meV, and 40 mev, respectively, (c) Photon 
energy vs. peak position in the chemical potential \i and (d) corresponding peak photocurrent at 
different temperatures. G = 40 meV is used for (c) and (d). 
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